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Abstract 



A model assuming coherent quadrupole-octupole vibrations and rotations is ap- 
plied to describe non-yrast energy sequences with alternating parity in several even- 
even nuclei from different regions, namely ^^^'^^^Sm, 154,156,158qj^ 236-jj ^,^^ i'^Oj^q^ 
Q \ Within the model scheme the yrast alternating-parity band is composed by the mem- 

5 ' bers of the ground-state band and the lowest negative-parity levels with odd angular 

momenta. The non-yrast alternating-parity sequences unite levels of /3-bands with 
higher negative-parity levels. The model description reproduces the structure of the 
►-j^ considered alternating-parity spectra together with the observed B(E1), B(E2) and 

cn ' B(E3) transition probabilities within and between the different level-sequences. B(E1) 

^^ , and B(E3) reduced probabilities for transitions connecting states with opposite par- 

ity in the non-yrast alternating-parity bands are predicted. The implemented study 



outlines the limits of the considered band-coupling scheme and provides estimations 
(3 ■ about the collective energy potential which governs the quadrupole-octupole proper- 

Cn ' ties of the considered nuclei. 



PACS: 21.60.Ev, 21.10.Re, 27.70.+q, 27.90.+b, 27.60. +j 



1 Introduction 

A typical manifestation of the reflection-asymmetric quadrupole-octupole deformation in 
the energy spectra of even-even atomic nuclei is the formation of level sequences with al- 
ternating parities [1]. Usually the levels with opposite parity are related through enhanced 
electric El and/or E3 transitions. The negative-parity sequence is shifted up with respect 
to the positive-parity sequence due to a tunneling of the system between the two opposite 
orientations along the principal symmetry axis. The magnitude of the energy shift corre- 
sponds to the softness of the shape with respect to the octupole deformation. The typical 
alternating-parity band is formed by the members of the ground-state (g) band and the 
levels of the lowest negative-parity sequence with odd angular momenta. In the relatively 
narrow region of the light actinide nuclei Rn, Ra and Th these two sequences merge into a 
single rotation band also called "octupole band" [21 |3l H] . The octupole band develops in 
the higher angular momenta and indicates the appearance of a quite stiff octupole deforma- 
tion. Away from the light actinide region both sequences diverge and do not form a single 
rotation band in the conventional meaning. Nevertheless, in some heavier actinides, like U 
and Pu and some rare-earth isotopes like Nd, Sm, Gd and Dy, they still remain related by 
El and E3 transitions, which indicates the presence of a soft octupole mode in the collective 
motion. In this case the term "alternating-parity band (or spectrum)" does not have the 
same strict meaning as in the light actinide nuclei but simply refers to sequences of levels 
with opposite parities which could be connected (coupled) through electric transitions. 

Various theoretical models have been developed over the years to explain and describe 
the formation of alternating-parity (or octupole) bands in the stiff and soft octupole regimes 
of coupling between the g-hand and the lowest negative-parity sequences in different nuclear 
regions [5]-|2n]. Particularly, a collective model assuming coherent quadrupole-octupole vi- 
brations and rotations ^ISJ was applied to the nuclei ^^°Nd, ^^^Sm, ^^^Gd and ^^^Dy with the 
presence of a soft octupole collectivity. Although the ^f-band and the lowest negative-parity 
bands in these nuclei were successfully described as members of an yrast alternating-parity 
band together with the attendant B(E1) and B(E2) transition probabilities, a question 
arises about the validity of such a consideration with respect to the higher-energy (non- 
yrast) part of the spectrum. 

The purpose of the present work is to clarify the above question within the model of 
Coherent Quadrupole-Octupole Motion (CQOM) |I8] by examining the possible formation 
of non-yrast alternating-parity structures in addition to the yrast band. For this reason 
the model scheme is extended by assuming that the excited /3-bands can be connected to 
higher negative-parity sequences with odd angular momenta. Therefore, it is supposed that 
the quadrupole-octupole structure of the spectrum develops along the non-yrast regions of 
the energy spectrum. Such a study provides not only a test of the model in the higher 
energy parts of the spectra, but also gives an interpretation of a larger number of data 
that may guide the experimental search for similar level structures in other nuclear regions. 
In principle, the systematic analysis of the non-yrast levels with alternating parity may 
favour different band-coupling schemes in the different nuclear regions allowing one to 
compare the capabilities of various theoretical models. For example, an extended study 
of non-yrast energy sequences with different parities has been implemented within the 



Extended Coherent States Model [16] by considering a coupling of the /3 and 7 bands with 
respective bands possessing the same spins but opposite parities, as well as a coupling 
between K'^ = 1"*" and K^ = 1~ energy sequences. In the model scheme of the present 
work the positive-parity /3-band appears connected to a negative-parity non-yrast sequence 
with odd angular momenta in the same way as in the yrast alternating-parity configuration. 
This is a consequence of the assumed mechanism of coupling between the quadrupole and 
octupole vibration modes. Therefore, the present work suggests a different band-coupling 
scheme and supposes a persistent role of the quadrupole-octupole motion in the forming 
of the higher-energy (non-yrast) part of the spectrum. Of course, by developing such an 
approach one should keep in mind the non- conventional meaning of the term "alternating- 
parity band" mentioned above. Also, presently the CQOM model is limited to excitations 
associated with the axial quadrupole and octupole degrees of freedom. Therefore, the study 
is focused on the related part of the collective spectrum, while other kinds of excitation 
modes as the 7-vibrations remain beyond the present consideration. 

The paper is organized as follows. In Sec. 2 the CQOM model is presented and the 
model mechanism for the appearance of non-yrast alternating parity bands is shown. Model 
expressions for reduced B(E1), B(E2) and B(E3) transitions in the non-yrast spectra are 
given in Sec. 3. In Sec. 4 numerical results and discussion on the application of the model 
to the nuclei of different regions are given. Sec. 5 contains concluding remarks. 

2 Model of Coherent Quadrupole— Octupole Motion 



The CQOM model [18] is a particular realization of the more general geometric concept of 
collective nuclear motion characterized by the quadrupole-octupole shape deformations ^ . 
The expansion of the surface radius R{9, (p), in polar coordinates, with respect to spherical 
harmonics up to multipolarity A = 3 is given by 



R{0, v) = Ro 



A=2 /x=-A 



where Rq is the spherical radius and oa^ are the twelve quadrupole and octupole collective 
coordinates in the laboratory frame. The collective coordinates are transformed into a 
body-fixed frame 

aA. = $^«AM^i.W, (2) 

determined by the "canonical" quadrupole coordinates Qq = 020 and 02 = 022 = 02-2 
and the three Euler angles 9 = (6*1, 6*2, 6*3). The remaining seven octupole coordinates a^^ 
(/i = —3, ...,3) together with oq and 02 determine the quadrupole-octupole shape of the 
nucleus. In the particular case of axial symmetry the quadrupole-octupole deformation 
represents a pear-like shape determined by the only non-zero coordinates P2 = 0,0 and 
^3 = 030. The respective physical states of the nucleus in the intrinsic (body-fixed) frame 
are characterized by the symmetrization group Dqo which consists of arbitrary (infinite 
number) rotations about the intrinsic z-axis and rotations about the axes perpendicular to 



z through the angle vr. In principal the symmetrization group of the nucleus in the intrinsic 
frame is determined by the rotations g satisfying a set of equations in the form D^^i^g) = 0, 



which in the case of axial symmetry is D]iQ{g) = for all yU 7^ [2T] . 

In the CQOM model [TB] the geometric concept is implemented in the limits of the 
axial symmetry. It is considered that the even-even nucleus can oscillate with respect to 
the quadrupole /52 and octupole P^ axial deformation variables, which are mixed through 
a centrifugal (rotation- vibration) interaction. The collective Hamiltonian of the nucleus is 
then taken in the form 

where 

W..A./) = ic.ft^ + ic,ft= + ;y||^. (4) 

with X{I) = [do + I {I + l)]/2. B2 and -B3 are effective quadrupole and octupole mass 
parameters and C2 and C3 are stiffness parameters for the respective oscillation modes. 
The quantity j7'(qiad+oct) _ (^d^f]^ _|_ d^f]^^ can be associated to the moment of inertia of an 
axially symmetric quadrupole-octupole deformed shape [23] with d2 and d^ being inertia 
parameters. The energy potential (jl]) represents a two-dimensional surface determined by 
the variables P2 and P3 with an angular-momentum-dependent repulsive core at zero defor- 
mation (see Fig. 1 in [18j). The parameter do in the centrifugal factor X{I) characterizes 
the repulsive core at / = and determines the overall energy scale for the rotation part of 
the energy. 

The model Hamiltonian ([3]) represents a D^o invariant. Also, it is important to remark 
that (|3]) corresponds to a class of axial-symmetric Hamiltonians [9] , [10] , [13] whose kinetic 
vibration parts are derived by ignoring the non-axial degrees of freedom (e.g. 7- vibrations) 
in a way similar to the approach of Davidov and Chaban [22] . The scalar product in 
the space of the wave functions (e.g. see Eqs. (2) and (4) in [13]) corresponding to the 
particular form of the (32- and /^s-derivatives in ([3]) is characterized by a unit weight factor, 

i.e. ($2|$i) = //t//32t;/33$^(/32,/33)$i(/32,/33). 

If a condition for the simultaneous presence of nonzero coordinates (/3™^'^, /3™°) of the 
potential minimum is imposed, the stiffness and inertial parameters are correlated as 
d2/C2 = d-i/Cs (see Eqs. (3)-(6) in [IS])- In this case the potential bottom represents 
an ellipse in the space of P2 and P^ which surrounds the infinite zero-deformation core (see 
Fig. 3 in [IE]). If prolate quadrupole deformations /92 > are considered, the system 
is characterized by oscillations between positive and negative /Ss-values along the ellipse 
surrounding the potential core. By introducing polar-type of curvilinear or, more precise, 
ellipsoidal variables 



V 
such that 



2{d2Pl + ds(3l 



d2 + ds 



and (h = arctan 




(32= prj cos (j), (3^ = qr] sin ((), (5) 



with 




V "2 

the potential (j4]) appears in the form 

Uiiv) 




and d 



l-^'^W 



\{d2 + d3 



(6) 



(7) 



where C = {d/d2)C2 = id/d3)Cs. 

Further, it is assumed that the quadrupole and octupole modes are represented in the 
collective motion with the same oscillation frequencies U2 = 003 = oo, with 





OJ 



The condition ([8]) imposes certain correlations between the mass, stiffness and inertia pa- 
rameters of the model Hamiltonian ([3]), corresponding to a coherent quadrupole-octupole 
motion of the system. Note that here the term "coherent" is used in the context of the 
mixing between the quadrupole and octupole degrees of freedom, which is different from 
the meaning of the same term used in [16]. In this case the Hamiltonian is obtained in a 
simple form 



Hnn = — TTTT 



<-qo 



2B 



d^ 
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dr]"^ T] dr] rf dcj)^ 



+ Uiiv) 



It allows an exact separation of variables in the wave function $(77, 0) 
subsequent equations for '(/'(r/) and (/?((/>) 

^' ,/ N 1 -9 ,, , 2B 



E- 




^iv) = 


-- 0; 






M<t^) - 


= 0, 



(9) 
ip{ri)ip{(f)) with the 

(10) 

fin 



where k is the separation quantum number. Eq. ( ITOl) with the potential ([7]) is similar to 
the equation for the Davidson potential [24j and has the following analytic solution for the 
energy spectrum [18] 

(12) 



En,k{I) = nw 2n + l + ^/WTbX(7) 



where u is defined in ([8]), n = 0, 1, 2, ... and b = 2B /{ti^d). The quantum numbers n and k 
have the meaning of "radial" and "angular" oscillation quantum numbers, respectively. The 
normalized "radial" eigenfunctions il){ri) are obtained in terms of the generalized Laguerre 
polynomials 



^Lkiv) 



2cr{n + 1) 

r(n + 2s + l] 



-^'''\c^yLl\cr^') 



(13) 



with c = VbC/H and s = {l/2)^/k^ + bX{I). Eq. ^ in the "angular" variable is 
solved under the boundary condition ip{—7i/2) = ip{7i/2) = 0. This is equivalent to the 
consideration of an infinite potential wall at /32 = (or ip = ±7r/2). Then one has two 
identical solutions for (32 > and /32 < 0. As mentioned above the physical space of the 
model is taken in the prolate /32 > half of the (/92,/93)-plane. (See Figs. 4 and 5 in [18] 
and the related text in that reference). Within this half-plane Eq. 01 ip has two different 
solutions with positive and negative parities, vr = (+) and n = (— ), respectively 



^^(0) = 


= 1/2/7? COs(fc0) , 


rC — J-, 0, O5 


¥>k{<P) = 


= ^/2/nsm{k(P) , 


A; = 2, 4, 6, 



(14) 
(15) 

Note that the square root term in the wave function ipikiv)^ Eq. f lT3|) . differs from the 
respective term used in Eq. (24) of Ref. |IS] by the factor c which is newly included in the 
numerator. (In \TE\ the quantity c is denoted by 'a' which in the case of odd nuclei leads 
to confusion with the notation for the decoupling factor.) One can easily check that this 
factor is necessary to normalize ipi kiv) to unity. The results for the transition probabilities 
obtained in [18] are not affected by the missing factor c due to the use of overall scaling 
factors in Eqs. (46) and (47) of [IS]. 

Since the consideration is restricted to axial deformations only, the projection K of the 
collective angular momentum on the principal symmetry axis is taken zero. Then the total 
wave function of the coherent quadrupole-octuple vibration and collective rotation of an 
even-even nucleus has the form 

KkiMoiv, 4>) = sj^^Di,,{e)^l,,{r^, 0) , (16) 

where DI,jq{6) is the Wigner function defined according to the phase convention in [25]. 
Note that due to a different phase convention in some other works, e.g. in [26] and |27j . 
the complex conjugated D-function appears in the rotation part. The relations between 
the different definitions of the D-function are given in Table 4.2 in [28]. The quadrupole- 
octupole vibration part of flTB]l is 

KkiivA) = i^LivMi^l^) ■ (17) 



The quantum numbers of the quadrupole-octupole vibration function fITTj) are determined 
by the requirement for a conservation of the T^P-symmetry of the total wave function 
( fT6l) . {V is the parity operator and TZ represents a rotation by an angle vr about an axis 
perpendicular to the intrinsic z-axis) The 7?.-symmetry of the rotation function D\jq{9) 
is characterized by the factor (—1)''^, while the action of V on (^'^j^jijq^cj)) gives the factor 
vr = ±. Then the conservation of the T^.'P-symmetry is equivalent to the conservation of 
the so called simplex quantum number simplex = 7r(— 1)^ = 1. This condition imposes a 
positive parity for the states with even angular momentum, and negative parity for the odd 
angular momentum states, i.e. one has 

Kkiiv, 0) = '^ikiv)vti(P) for / = even 
Kkiiv, 0) = ^ikivWki(f>) for I = odd. 
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It should be noted that the above conditions are in conjunction with the transformation 
properties of the variables r] and in (fT7|l under the rotation TZ (77 is invariant, while 
changes in sign) so that together with the simplex conservation condition the total wave 
function (ITB]1 appears to be an Dqo invariant as it should be due to the axial symmetry. 

The structure of the energy spectrum is determined by the oscillator quantum numbers 
n ("radial") and k ("angular") in Eq. (IT^ . Since according to Eqs. (Tn|) and (TTB1) k obtains 
different values for the states with opposite parity the energy sequences with even and odd 
angular momenta corresponding to a given n appear shifted to each other, i.e. a parity shift 
effect is observed. In \XS] it was supposed that the ^f-band and the lowest negative-parity 
band belong to a n = set with k = k^^' = 1 ioi g and k = k^' = 2 for the negative- 
parity band. In the present work the model scheme is extended through the following three 
suppositions. 

i) The energy spectrum determined by the coherent axial quadrupole-octupole vibrations 
and rotations consists of couples of level-sequences with opposite parity. The sequences in 
each couple are characterized by the same value of the quantum number n = 0, 1, 2, ... and 
by different values oi k, k = kn = lor 3 or 5 or ... for the even- J sequence, and k = kn = 
2 or 4 or 6 or... for the odd-/ sequence. 

ii) The lowest values of the "radial" quantum number n correspond to the lowest alter- 
nating parity bands, with n = being the yrast band, n = 1 corresponding to the next 
non-yrast alternating parity structure and so on. The values of the "angular" quantum 
number k are not restricted and should only satisfy the parity condition in i). The par- 
ticular values of ki and kn can be determined so as to reproduce the experimentally 
observed parity shift in the set of levels with a given n. 

iii) Due to the coherent interplay between the P2 and P^ variables in the oscillation mo- 
tion, the excited /3-bands in even-even nuclei can be interpreted as the positive-parity coun- 
terparts of higher negative-parity sequences, or as the members of non-yrast alternating- 
parity bands. 

Based on the above assumptions the extended alternating-parity spectrum of an even- 
even nucleus can be considered in the following form. 

Yrast alternating-parity set (n = 0): unites the ^f-band {k = /cq ) /J = 0f,2f,A^,Q^,... 
with the first negative-parity band denoted here as nl {k = k^ ) /^ = lf,3f,5]',...; 

First non-yrast set (n = 1): unites the first /3-band denoted by h\ {k = k\ ) I'^ = 
J, 2 J, 4^, ... with the second negative-parity band denoted hy n2{k = k\ ) /^ = 1^, 3^, 5 
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S+)\ T7T _ n+ 



Second non-yrast set (n = 2): unites the second /3-band 62 {k = k^ ) /J = 03,23,43,... 
with the third negative-parity band ra3 [k = ^2 ) K — l-T' 3J, 5J, ..., and so on, higher 
non-yrast sequences, where v = 1,2,3,... is the consequent number of the appearance of 
a state with a given angular momentum. Also, it is convenient to use the band labels 
introduced above to denote the different excited states as for example 2+, 1"^^, 0^^, 1~2 etc. 
Obviously the above model scheme makes no claim to exhaust the entire collective spec- 
trum but rather provides a tool to identify the extent to which the considered quadrupole- 
octupole motion can influence the excited band structures in even-even nuclei. In the end 



of this section, it should be remarked that the extension of the model to higher energy 
levels, together with assumption ii), which releases k from the fixed values k^~^^ = 1 and 
k^~^ = 2 (originally imposed in [18] for the yrast case), now requires a new readjustment of 
the model parameters. 

3 Transition probabilities in the non- yrast quadrupole- 
octupole states 

As the B(E1) and B(E3) reduced transition probabilities are known to provide a sensitive 
test for the structure of the alternating-parity sequences it is of special importance to 
examine their behaviour in the non- yrast part of the spectrum. The basic CQOM concept 
for the electromagnetic transitions has been given in [18] . Here the formalism is modified so 
as to describe B(E1), B(E2) and B(E3) reduced transition probabilities in the higher lying 
alternating-parity bands along with the extended treatment of the model energy quantum 
numbers. A more essential modification is related to a generalization of the angular part 
of the electric transition operators dictated by the complicated quadrupole-octupole shape 
density distribution inherent for the coherent motion mode (see below). In addition the 
E1-E3 charge factors are treated explicitly and the model parameters p and q (jH]) providing 
information about the potential shape are considered without including them into scaling 
constants. 

The reduced transition probability for an electric transition with a given multipolarity 
A between model states ( !T6l) with n = rii, k = ki, I = Ii and n = Uf, k = kf, I = If is 

The operators for electric El, E2 and E3 transitions have the following general form 



2A -I- 1 
M^{EX) = .I ^^ ^^ , Qao^o\, A = 1,2,3, /i = 0, ±1, ..., ±A. (19) 



47r(4 - 36, 



A,l 



The vibration parts of these operators are given up to the first order of (32 and (S^, for E2 
and E3, and in second order, for El, as 

Qio = Ml/32/33 (20) 

Qxo = Ma/3a, a = 2, 3. (21) 

The electric charge factors Mx (A = 2, 3) are taken as [29] 

^^ = ^T^^^T^^^^^o ' ^ = 2, 3 , (22) 

a/(2A+ IJTT 

where Rq = tqA^^^, rg ~ 1.2 fm, Z is the proton number, and e is the electric charge of the 
proton. The charge factor Mi is taken according to the droplet model concept [30]- [32] in 
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the form [TO 



9AZe^ /I 15 \ 
Ml = ^ - + r , 23 

where the quantities J and Q are related to the volume and surface symmetry energy, 
respectively and their values are assumed in the limits 25 < J < 44 MeV and ^7 < Q < 70 
MeV [33] (see also the values below Eq. (79) in [ID]). In the present work fixed average values 
of these quantities J = 35 MeV and Q = 45 MeV are used for all considered nuclei. One 
should remark that so far there is no unique approach to estimate the factor Mi. Therefore, 
here in fl23|) the proton charge e is replaced by an effective charge eLj- which is considered 
as an adjustable parameter and can be different from one. Note that to obtain the B(E1) 
transition probabilities in the units e^ -fm^, and subsequently in Weisskopf units one has to 
take into account that e^ = 1.4399764 MeV ■ fm [or eVMeV^ = (1.4399764)2e2 ■ fm^] which 
leads to an additional multiplication factor 1.4399764 in f l23|) when numerical values are 
produced. 

In the space of the ellipsoidal coordinates ([5]) , ([6]) one has 

Qio = Mipqrf cos (j) sin (f) (24) 

Q20 = M2pri cos (f) (25) 

Q30 = M^qr] sin (f). (26) 

The definitions of the operators fl2Ul) - fl?Il) and fl2^ - fl2Bl) originally correspond to a situation 
in which the nuclear shape is characterized by fixed values of the deformation parameters 
/32 and /Ss. In this case the density distribution of the collective state is characterized by 
a single maximum in the space of P2 and Ps. In the case of the model potential (jlj) taken 
with an elliptic bottom the density distribution can be characterized by more than one 
maximum. Indeed, the density of the model state fll7p is characterized by a different number 
of maxima depending on the quantum number k. This feature is a result of the assumed 
soft quadrupole-octupole mode. It is illustrated graphically in Appendix A, where the 
density distribution of the state flT7|) in the space of the quadrupole-octupole deformations 
is plotted for different fc-values at ra = after transforming the wave function ^^^j- in the 
(/32,/33) variables. It is seen that for /32 > the number of maxima is equal to k and 
by analogy with the acoustics may be interpreted as the number of "overtones" which 
characterize the coherent collective oscillations of the system. Thus, it appears that the 
transition operators should connect states with different numbers of maxima (or overtones). 
In the space of ellipsoidal variables the positions of the maxima are determined by the 
angular variable 0. On the other hand the original operators fl24|)- fl26j) do not take into 
account the presence of multiple maxima in the shape density distributions of the different 
states. One particular effect due to this circumstance is that the integrals over the angular 
part of fl26|) . sin 0, vanish when the difference between the k numbers of the initial and final 
states is larger than a unit and the respective B(E3) transition probabilities vanish too. 
This limitation is removed if the operators are generalized appropriately. The most general 
forms of the angular parts of the operators corresponding to the first orders of P2 and Ps 



according to (!25|) and (l26l) can be sought in terms of a Fourier expansion with respect to 
through the replacements 



COS0 — 7- v42o(0) = y^o-20 cos(/c0), siu — )■ ^30(0) = y^cLso sm{k(f)), (27) 



4? cos(/c0), sin0 -^ ^30 (0) = X] ''^''^ ' 
fc=i fc=i 



where the expansion coefficients should be chosen so as to provide an appropriate conver- 
gence. A choice made here for both type of coefficients is a^''^ = 1/k for which the above 
expansions can be obtained in analytic form 

A2o{4>) = f;^^4^ = -^[ln2 + ln(l-cos0)] (28) 



fc=i 

oo 



AM - E^-^ + 'Fl„o.(^ , (29) 

k=l ^ ^ 

where the Floor function maps a real number to the largest previous integer. Then the 
angular part of the second order operator (1241) can be generalized in an obvious way 

i • 1 A r i\ A f i\ A r i\ v~^ v^ cos(m0) sin(n0) .^^, 

cos0sm0 ^ Aio(0) = A2o(0)^o(0) = V V —^ — -^^- (30) 

^ — ^ ^ — ^ m n 

m=l n=l 

Note that the first terms of the above expansions represent the original angular parts in 
(l2^ -(l26l). So, the new angular operators (I28l) -( l30l) . which are extensions of the old ones, 
provide a connection between states whose "dynamical" deformations (i.e. the probabil- 
ity distribution in the deformation space) are characterized by the co-existence of a large 
number of maxima. These specific shape properties of the system are due to the assumed 
coupling between quadrupole and octupole degrees of freedom. Now the operators f l2^ -f l26|) 
are redefined as 

Qio(r/,0) = MmV^AM) (31) 

Q2o{r]A) = M2PvA2oicf>) (32) 

QsoivA) = MsqvAM). (33) 

After carrying out the integration over the rotation part in (llSp one obtains 

2A -I- 1 

B{E\-n,hh -^ Ufkjlj) = -— M^{)\{]\If{])^Rl{n,hh -^ nfylj), (34) 

47r(4 - 6dx,i) 

which involves the squares of the Clebsch-Gordan coefficient and the matrix element of 
the electric multipole operators (13T]) - (133|) between the quadrupole-octupole vibration wave 
functions (TT7|) 



Rxin^hh ^ njkjij) = (KUfiM<l^)\Qxo\K]kjSv,<f>)) ■ (35) 
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By further separating the integrations over the "radial" variable rj and the "angular" vari- 
able (j) in (135|) according to (TT7|) one obtains 

Riiriikili^ njkflf) = MipqS2{ni,Ii;nf,If)l'^"''^'{ki,kf) (36) 

R2{nikili -^ Ufkflf) = M2pSi{ni,Ii;nf,If)I^"'^^{ki,kf) (37) 

Rsiuikili -^ rifkflf) = MsqSiiui, If, Uf, Ij)!^"''^ {k,,, kf), 



(3^ 



where 



Si{ni,Ii;nf,If) = 


I 

POO 


S2{ni,Ii;nf,If) = 


f 



dvi^'4iv)vVn\iv) (39) 

dr]ij'4.{r])7]'^i^^{7]), (40) 



and 



C^'ih^kf) = ^l^[A,oi<P)^7^{cf>)^imd<P, A = 1,2,3. (41) 

The integrals over rj, (l39l) and (l40l) . involve the "radial" wave functions ( iT3l) . Analytic 
expressions for these integrals are given in Appendix B. The integrals over ( 14T1) involve 
the "angular" wave functions flT^ and/or flT5|) . The explicit forms of these integrals with 
the relevant parities ttj and ttj are given in Appendix C 

From the generalized definitions (!3T|) - (!33l) of the operators Qao it is seen that the inertial 
factors p, q and their product pq defined through Eq. (jS]) are not included in any scaling 
factors, as done in Ref. [IS] and can be considered as model parameters. Actually, p and q 
are not independent. From ([6]) it can be easily seen that 1/p^ + 1/g^ = 2. Then q can be 
expressed by p as 

q = , ^ with p> ^^ 0.7071. (42) 

y2F^I V2 

The inequality in ( |42|) corresponds to the condition (i2 < 2(i. Analogically one can express 
p hj q with the condition (^3 < 2d. (Note that for p = 1 one has q = 1 and pq = 1.) Here 
the adjustable parameter is chosen to be p. It should be noted that with the involvement of 
the new parameter p the scaling factors in Eqs. (46) and (47) of Ref. [IB] are not considered 
anymore and the charge factors M2 and M3 are directly calculated in fl22|) . Also the charge 
factor Ml is directly calculated in f l23|) . but with the effective charge elj-j being adjusted 
to determine the correct scale of the B(E1) transition probabilities with respect to B(E2). 
From another side the parameter p determines the relative scale between B(E2) and B(E3). 
It is interesting to remark that p does not play any role if the model energy levels are 
fitted without taking into account transition probabilities. However, in this case there is an 
ambiguity in the choice of the inertial parameters d2 and d^. This is seen by the following 
relations between the parameters of the potential (j4]) and the fitting parameters u and b, 
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imposed by the coherent condition ([8]) 

C, = ^d,, Cs = ^d, (43) 

B2 = \d2, 53 = ^4, (44) 

in which ^2 and (is are not determined. (The parameter d^ does not enter these relations.) 
This means that for a given set of parameters cj, h and d^ the energy spectrum corresponds to 
an infinite number of potential shapes with different eccentricities of the ellipsoidal bottom. 
Now, after determining the parameters p and c with respect to transition data one gets 

c^2 = ^, d^ = {2p^-l) — , with d=^. (45) 

p^ p-^ UJO 

Thus for given values of the parameters u, b, c and p the original parameters B2, -B3, C2, 
C3, d2 and d^ of potential (jl]) are fixed, and given additionally do, its form is unambiguously 
determined. 



4 Numerical results and discussion 

The extended CQOM formalism was applied to several nuclei, namely ^^^'^^'^Sm, i54,i56,i58q^^ 
^^^U and ^°°Mo, in which one or two non-yrast alternating-parity bands can be constructed 
by the experimentally observed /3 and higher-lying negative-parity levels. In these nuclei 
a number of data on El and/or E3 transitions are available providing the possibility to 
test the complete model scheme. In all selected nuclei the experimental data [31] provide 
well determined yrast and first non-yrast alternating-parity bands except for ^°°Mo where 
the structure of the non-yrast band is proposed here on the basis of the model analysis 
(see below). In three of the nuclei, ^^^Sm, ^^"^Gd and ^^^Gd second excited (non-yrast) 
alternating-parity bands are additionally considered. The structure of these bands is not 
clearly determined in the experimental data. Therefore, the model description and predic- 
tion provides a possible interpretation of the respective experimental levels. In this meaning 
the present description not only provides a test for the CQOM model scheme, but also sug- 
gests a possible classification of some highly non-yrast excited states whose interpretation 
in the experimental data bases is not unambiguous. 

The model description is obtained by taking the theoretical energy levels En,k{I) = 
Enkil) — E ,(+)(X)) from Eq. (fT2!l . The parameters w, 6, d^, c, p and the effective charge 

e\ff have been adjusted by simultaneously taking into account experimental data on the 
energy bands and the available B(E1)-B(E3) transition probabilities. The parameter values 
obtained in the considered nuclei are given in Table 1. The resulting values of the original 
Hamiltonian parameters in (|3]) and dl]) are given in Table 2. For each nucleus the calcu- 
lations are performed in a net over the values of the "angular" quantum numbers k with 
appropriate parity in the limits 1 < fc < 20. In all nuclei sets of values for the /c-quantum 
numbers providing the best model description of both, energies and reduced transition 
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probabilities, are obtained. These values are given in Figures 1-7 where the theoretical 
and experimental energy levels of the considered nuclei are compared. The theoretical and 
experimental values of the B(E1), B(E2) and B(E3) transition probabilities are compared 
in Table 3. Model predictions for some not yet observed transitions are also given there. 

The results in Figs. 1-7 show that the model scheme correctly reproduces the structure 
of the alternating-parity spectra in the considered nuclei with a reasonably good agreement 
between the theoretical and experimental energy levels. The correct reproduction of the 
mutual displacement of the different positive and negative parity sequence is related to 
the involvement of k quantum number values larger than 1 and 2. On the other hand 
the determination of the fc-values is strongly dictated by the interband transitions between 
the positive- and negative- parity levels as well as by the transitions between the different 
alternating-parity sequences. The above remarks explain why for the nuclei ^^^Sm and 
^^^Gd new sets of /c-quantum numbers appear together with renormalized values of the 
fitting parameters compared to the previous descriptions limited to the yrast bands [18] 
(see below). One should remark that at the same time the main (radial) oscillator quantum 
number n is uniquely determined, n = for the yrast sequence, n = 1 for the first excited 
alternating parity band and so on, as explained in the end of Sec. 2. 

In ^^^Sm the yrast band is described together with the first excited band (see Fig. 1). The 
calculations provide two identical couples of k values {k^~^^ = 1, k^~'^ = 8) for each band. 
Thus it is seen that fc*^~^ obtains a value larger than the lowest even value 2 considered 
in [TB]. From Table 3 one can see that with this configuration of fc- numbers the model 
fairly good reproduces the data [35] on the B(E2) intraband transition probabilities in the 
ground-state band (g) and on the B(E1) probabilities for transitions between the g- and 
the first negative-parity band (nl). Some interband E2 transitions between the g- and 
the first /9-band (61), as 2'^^ — )• 2+ and 4^^ — )■ 4+ are also well described, while others like 
4^^ — !■ 6J" are overestimated. The calculated intraband transitions in the 61-band are in 
rough agreement with the experimental data, while the El intraband transition 1"^^ — )■ 2^^ 
is overestimated by an order. The E3 transition probability B{E3; 3^^^ ^0^) = 14 W.u. 
[36] is exactly reproduced due to the adjustable parameter p which determines the factor q 
in fl38l) according to Eq. (H2|) . This allows one to predict other E3 transitions like l„i — > 4^, 
and other similar transitions between the 61-band and the second negative-parity band 
(n2) as shown in Table 3. Although not all theoretical transition probabilities are in strict 
agreement with the experimental data it is seen that the model scheme correctly takes into 
account the different scales of the various kinds of probabilities. A similar behaviour of 
transition probabilities is observed in the other considered nuclei. 

In ^^^Sm totally three alternating-parity bands are considered as seen from Fig. 2. The 
positive-parity states of the second excited band are interpreted in [3^ as members of a 
second K'^ = 0"*" band, or of a second /3-band (62). The respective negative-parity levels are 
selected in the present work among levels for which there is no interpretation given in [31] . 
Here they form a third negative-parity band (ra3). From Table 3 it is seen that the intraband 
B(E2) transition probabilities in the g-hand of this nucleus are reasonably well described 
up to / = 10, while the B{E2] 12+ — )■ 10^) value is considerably overestimated. The B(E1) 
probabilities between the nl- and g-hands are also well described. The theoretical interband 
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transition value B(£'2;0^ — ;■ 2+) is by an order smaller than the experimental one. For 
the other similar E2 transitions like 2^^^ — )■ OJ" the theoretical values are obtained below the 
upper limits given for the respective experimental data [35] . 

In ^^^Gd, again, three alternating-parity bands are considered. The model description 
is given in Fig. 3. Here, the second non-yrast band is constructed by the second excited 
K'^ = 0^ band and a 3~ state with energy 1796.96 keV ^34j- Although the latter is in- 
terpreted in [3lj as a member of a K'^ = 2~ octupole band it reasonably fits the present 
scheme as a member of an n3-sequence. In this nucleus the B(E1)-B(E3) transition prob- 
abilities are also reasonably described with the largest discrepancies between the theory 
and experiment, about a factor of 2, being observed for the E2 transitions 2^^ — )■ 0^^^ and 
0^^ — ;• 2J" (see Table 3). Note that here the theoretical B(E1) value for the interband tran- 
sition B(£^l; 1~^ — )■ 2^^) = 0.0064 W.u. is obtained close to the experimental one, 0.0099 
W.u. 

In ^^^Gd two alternating-parity bands, the yrast and first excited, are considered (see 
Fig. 4). The B(E2) and B(E1) transition probabilities between the members of the g-, bl- 
and nl-bands are well described with a few exceptions as in the transitions 4^^ — )■ 2^ and 
4^]^ — )■ 2i" for which the B(E2)- values are underestimated with respect to the experiment by 
a factor of about two and an order, respectively (see Table 3). On the other hand the model 
predictions for the B(E1) transition probabilities between the second negative-parity band 
n2 and the f^-band suggest 2-3 orders of magnitude in suppression compared to experimental 
data. 

In ^^'^Gd, three alternating-parity bands are considered (see Fig. 5). Similarly to ^^'^Gd 
the 1~ and 3" states included in the second excited band enter the present model scheme 
as n3- members, while in [31] they are interpreted as members of a K^ = 1^ octupole 
band. For this nucleus, quite a large number of data on B(E1) and B(E2) transition 
probabilities are available [55]. One should remark that compared to the other considered 
rare-earth nuclei ^^^Gd is closer situated to the region of pronounced rotation collectivity. 
From Table 3 it is seen that the theoretical intraband B(E2) probabilities in the gf-band 
of ^^®Gd faster increase with the angular momentum compared to the experimental data. 
On the other hand six experimental B(E1) values for the transitions between the g- and 
the 77.1-bands are described quite well. It is remarkable that an experimental estimation 
for a El transition between the n2- and 61-band is available with B(ii^l;3~2 — ^ 2^^) > 
0.00035 W.u. This circumstance is in a conjunction with the model assumption about 
the quadrupole-octupole coupling of both bands. The model description predicts for this 
probability a smaller value of 0.00011 W.u., which is of the same order as the B(E1) 
values connecting the g- and nl-bands. Further, model prediction values for similar B(E1) 
transition probabilities as B{E1; l'^ -^ 0+ ) = 8 x 10"^ W.u. and B{E1; 1'^ -^ 2+ ) = 0.0002 
W.u. are given in Table 3. Also there one can find available experimental estimations for 
intraband transition probabilities like B(i?2; 5^^^ — )■ 3^^) = 369 W.u. and B(i?2; 3^2 ~^ 
1,~2) > 1600 W.u. which are underestimated by the theory. In addition, a number of B(E1) 
transition probabilities from 61- to nl-, from n2- to g- and from 62- to nl- and n2-bands 
are generally underestimated by one or two orders of magnitude. 

In ^^^U two alternating-parity bands, the yrast and first excited, are considered (see 
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Fig. 6). This nucleus was selected because of the possibility to examine two observed 
reduced probabilities for E3 transitions, namely B{E3; 1~^ — )■ 4^) = 62 W.u. [35] and 
B(£^3;3„i — i- Og) = 22.9 W.u. [36j. From Table 3 it is seen that the first one is exactly 
reproduced. The second one is underestimated by the theoretical value, 15 W.u., which is 
still reasonably close to the experiment. For the B(E2) transition probabilities within the 
g-hand the description is good as overall up to a quite high angular momentum J = 26. The 
experimental B(E1) value for the transition 1^^ — )■ 0^ is exactly reproduced because of the 
use of the effective charge. This allows one to predict other B(E1) transition probabilities 
in the described spectrum which are given in Table 3. 

In ^*^°Mo the experimentally observed 2^ state with energy 1463.9 keV is considered in 
[M] as a possible member of a /3-band. However, the present scheme suggests that the 2"'" 
state belonging to this band should lie essentially lower. The calculations show that the 
experimental 2^ state with energy 1063.78 keV considered in [3lj as a possible member of 
a 7-band is more appropriate as a /3-band member. The result in Fig. 7 shows that if this 
state is included in the 61-band (in the present notations) a non-yrast alternating-parity 
sequence can be constructed and reasonably well described by taking three additional states, 
namely 1~ at 2156 keV, 3~ at 2369.6 keV and 4"*" at 1771.4 keV from the set of available 
but not interpreted data for ^"^^Mo [31]. The observed B(E1)-B(E3) transition probabilities 
are reasonably described as seen from Table 3. The main discrepancy between the theory 
and the experiment, a factor of 5, is obtained for the E2 intraband transition 2^^ — )■ O'l^. 

The following comments on the model results can be made here. The parameters of the 
fits shown in Table 1 reflect the common collective structure of the various energy sequences 
{g, 61, 62, nl, n2, n3) in a given nucleus, while the sets of k values given in Figs 1-7 reflect 
their mutual dispositions. Note that the parameters for ^^^Sm and ^^^Gd are essentially 
renormalized compared to the fits of the yrast band only [TB]. As seen from Table 1 
the parameters u and b, which are responsible for the rotation-vibration behaviour of the 
different sequences, vary relatively smoothly between the different nuclei. The parameter 
do, which is responsible for the shape of the potential at zero angular momentum, shows 
more pronounced differences in its values, especially for the nuclei from different regions as 
^^^U and ^°°Mo. Also, the parameter c, which determines the overall scale for the transition 
probabilities in the "radial" integrals, considerably varies, while the parameter p which is 
related to the quadrupole and octupole contributions to the moment of inertia changes 
quite smoothly. It is remarkable that in three nuclei, ^^^Sm, ^^^Sm and ^^^Gd, the effective 
charge for the El transitions is practically unit which means that there is no need for this 
parameter to describe them. In ^^^Gd it is still close to 1, while in the other three nuclei 
its need for the model description is already essential. 

By using the relations (IH|) and (1451) between the model parameters in ellipsoidal coor- 
dinates and the parameters of the original Hamiltonian, ([3]) with (jlj), one can obtain the 
latters from the values given in Table 1. Subsequently one can obtain the semi- axes (sa) 
(3^^ and (3f^ of the eUipsoidal potential bottom in the space of the quadrupole-octupole 
variables given by 

/3r(/) = [2XiI)/d,C,]'/\ A = 2, 3. (46) 

(For more details see the text after Eqs. (3) and (4) of [H].) The resulting values of the 
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parameters B2, B^, C2, C3, d2, d^ and the semi-axes are given in Table 2. Note that in the 
present work they are not directly adjusted, but obtained as a result of the adjustment of 
the parameters w, 6, rfo, c, p, and e\fj. As such they only give a rough estimation about 
the order of the potential parameters and its shape. One can see that for ^^^'^^^Sm and 
i54,i56Qj-| w^q^q parameters vary relatively smoothly, while for the remaining three nuclei 
they show some essential fluctuations. The values of the P^ semi-axis are obtained close 
to the known values of the static quadrupole deformations in these nuclei while the values 
of the octupole semi-axis P^ appear considerably larger. This result is correlated with 
the larger values of the quadrupole stiffness parameters C2 compared to the values of C3. 
Hence the present parameters correspond to a vibration motion with a larger softness of 
the system with respect to the octupole mode compared to the quadrupole one. A closer 
look on the formalism shows that the ratio between both semi-axes is related to the matrix 
elements of the quadrupole and octupole electric multipole operators (132|) and (133|) . By 
using (03]), (USD and (USD in (gnD one finds that 

^ = ? = ^— . (47) 

It is seen that the ratio P^/ Pf" depends on the inertia factors p and g, Eq ([6]), which 
determine the strength of the E2 and E3 transitions, respectively. This ratio is less than 
1 for p > 1 (g < 1). It can be easily checked that to obtain (if'/P^ < 1 one has to 
introduce an additional scaling constant, C3, having the meaning of an effective charge for 
the octupole mode. Then the octupole charge factor is renormalized as M3 = C3M3. The 
numerical analysis shows that if C3 is chosen in the limits 2 < C3 < 4 the parameter p is 
renormalized so that q — )■ g/3 and the same theoretical levels and transition probabilities are 
obtained with 13^ < 13 f in correspondence to the usually observed values of the deformation 
parameters P2 and P3. For example if C3 = 4 one obtains the following set of renormalized 
parameters for ^^^Gd, c' = 269.6, p' = 1.197, eLr = 1.512, while the parameters u, b and 
do remain unchanged compared to the values given in Table 1. Compared to the values in 
Table 2 the renormahzed parameters for ^^^Gd are B3 = 1146 ^V^eV, C3' = 108 MeV, 
ds = 777 ^^/MeV and Pl'^' = 0.192, while the other parameters referring to the quadrupole 
deformation remain unchanged. It is seen that now the length of the potential bottom semi- 
axis in the /33-direction corresponds to a more realistic octupole deformation. This result is 
equivalent to the involvement of a renormalized octupole operator Q'^Qir], 0) = csQso^r], 0). 
Since the use of such an effective charge does not change the model description but only 
leads to the renormalization of the parameters it is not considered in the present work. 

Further, it is important to comment the obtained configurations of quantum numbers 
kn and kn which characterize the energy shifts in the described alternating-parity spec- 
tra. From Figs. 1-7 it is seen that the relevant energy shift in the excited level sequences is 
obtained through a jump of k over several lower values. In this way certain low-lying states 
available in the scheme do not enter the considered spectrum, while others lying at higher 
energy are used to obtain the model description. This result is a consequence of the fact 
that the same oscillation frequency u is imposed to all alternating-parity bands. Actually, 
the non-yrast bandheads and the energy shifts could be reproduced through the lowest pos- 
sible fc- configurations [kn = 1, kn = 2] ii separate vibration frequencies are considered in 
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the different bands. Speaking about A; as a number of angular oscillation quanta (phonons) 
it appears that the restricted freedom of the frequency imposed by the coherent condition 
is compensated in the model description by the presence of a larger number of quanta on 
which the rotation bands are built. Since the eventual consideration of different oscillation 
frequencies would correspond to the introduction of parameters external for the model the 
larger numbers of quanta are retained in the present work. The obtained pairs of values 
kii and kn for the quantum number k provide a detailed systematic information about 
the mutual disposition of the positive- and negative-parity bands in the different nuclei, 
and subsequently, about the evolution of the quadrupole-octupole spectra in a given nu- 
clear region. It should be noted that the involvement of the extended transition operators 
(l3Tl)-(l33l) in the present CQOM development is related to the appearance of larger fc- values 
and the subsequent large ^-differences taken into account in the electric transition proba- 
bilities. These features of the model can change if it is applied beyond the coherent-mode 
assumption. In this case the unrestricted Hamiltonian ([3]) can be diagonalized by using the 
present analytic solution as a basis. Then the parameters in ([3]) can be directly adjusted 
to describe the spectrum without restriction of the quadrupole and octupole oscillator fre- 
quencies. This could allow one to construct the spectrum by always choosing the lowest 
possible eigenvalues, while the structure of the spectrum obtained in the present analytic 
solution could only guide the construction of non-yrast bands. Work in this direction is in 
progress. 

Finally, it should be noted that the present model descriptions are obtained within some 
natural limits of the applied formalism with respect to experimental data. It is well known 
that rotation terms like the one entering the model potential can only describe smooth 
changes of the rotation spectra with increasing angular momentum, as for example the so 
called "centrifugal stretching" . The treatment of angular momentum regions where sharper 
changes in the rotation spectrum due to changes in the intrinsic structure like backbending 
effects occur, needs a special development which is not the subject of the present work. That 
is why in some of the considered nuclei descriptions and/or predictions of rotation levels 
with very high angular momenta are avoided, especially in the cases where the negative- 
parity levels are not observed. An exception is done for ^se-jj {Yig. 6), where higher-spin 
negative-parity levels were predicted in accordance to the last observed state with even 
angular momentum. This prediction should be meaningful since in the actinide region the 
rotation spectra exhibit more regular rotation motion in the high-spin regions. On the 
other hand the prediction of missing low-spin states, like the 1~3 level in ^^^Gd and the 6^^^ 
and 5~2 levels in ^"^^Mo, as well as, a number of not observed transition probabilities shown 
in Table 3 should be also reasonable in the present framework. In this meaning the applied 
CQOM model scheme rather describes the "horizontal" evolution of the alternating-parity 
spectra beyond the yrast line than the high-spin properties of individual rotation bands. 

5 Concluding remarks 

The present work provides a model description and respective classification of the yrast and 
non-yrast alternating-parity spectra and the attendant B(E1), B(E2) and B(E3) transition 
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probabilities in several rare-earth nuclei, one U and one Mo nucleus within the collective 
model of Coherent Quadrupole and Octupole Motion (CQOM). The theoretical formalism 
and the obtained model descriptions outline a possible way for the development of nu- 
clear alternating-parity spectra toward the highly non-yrast region of collective excitations. 
In the considered scheme the different negative parity level-sequences appear in couples 
together with the ground-state band and the excited /3-bands. On this basis the model 
predicts possible E(l) and E(3) transitions between states with opposite parity within var- 
ious alternating-parity bands. The presence of experimentally observed E(l) transitions 
between such states in the non-yrast part of the spectrum is noticed. Further experimental 
measurements of electric transition probabilities would be very useful to check the possible 
coupling of non-yrast energy sequences with opposite parities. It was demonstrated that 
the considered scheme can be used for the interpretation of data on excitation energies 
whose place in the structure of the collective spectrum has not yet been determined. The 
approach was applied to selected nuclei for which a relatively large number of data on 
B(E1)-B(E3) transitional probabilities are available, but it can be easily extended to wider 
ranges of nuclei especially in the rare-earth and actinide regions. Further, the formalism 
takes into account the complex-shape effects in the motion of the system and in addition 
provides estimations about the shape of the quadrupole-octupole potential which governs 
the collective properties of the considered nuclei. More refined model descriptions and real- 
istic estimations about the potential shape can be obtained beyond the limits of the present 
coherent-mode assumption. Work in this direction is in progress. 
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Appendix A: CQOM shape-density distributions 

The density distribution of the CQOM vibration state in the space of the quadrupole- 
octupole shapes is given by the square of the wave function (TT71) . Pnfe/(/32, /Ss) = |'^nA:/(/^2, /Ss)^, 
after a transformation from the ellipsoidal coordinates (77, 0) to the deformation coordinates 
{P2, Ps)- In Fig. 8 three-dimensional plots of pnki are given for the lowest k = 1 and k = 2 
states for n = and for the schematic parameter values u = 0.3MeV/h, b = 3h~^, 
do = 100 h"^, 62 = 300 h'^/MeV, 4 = 500 h'^/MeV. Note that according to the discussion in 
the end of Sec. 3 the shape of the potential is determined unambiguously when the values 
of the inertia parameters d2 and d^ are given. In Fig. 9 two-dimensional plots showing 
the maxima of pnki for k = 1 — 4 are given together with contours showing the ellipsoidal 
potential bottom for the above set of schematic parameters. 
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Appendix B: Explicit form of the integrals over 77 

The integrals over t], ( l39l) and ( l40l) . can be written in the following common form after 
taking into account the explicit expression for the radial wave functions ( IT3l) 



Jo 

/"OO 

Jo 



(4J 



where 1 = 1,2, Si = {l/2)^kf + hX{Ii), Sf = {1/2) J k) + hX{If) and 



N = Nn„nf{c,Si,Sf) = 2C 



r{nf + l)r{ni + l] 



r(n/ + 2sf + l)T{n, + 2si + 1) 



(49) 



To derive an explicit expression for the integral ( HHl) one can apply the substitution erf = x 
with dx = 2cridri, such that 

^'^'^^ = ^^^'^'^^- (50) 



Then Eq. (H8l) reads as 

Nn„nf{c,Si,Sf) 



S,{n,,U;nf,If) = "'^J^/r' [e-^x-^^^^'^Ll'l{x)Ll-{x)dx. (51) 

By using known formulas for integration of two generalized Laguerre polynomials with 
different real ranks |37], [38] one obtains (|5T]) in the following explicit form 



Si{ni,Ii;nf,If) 

Nn,,nfic, Si, Sf) T{nf + 2s f + 1) Tjnj + gj - g/ - |) Tjsj + gj + | + 1) 
2ci+'/2 r(l + 2s/) r(si-s/-l) rzjn/! 

/ / / ^ . 

X 3F2 -ri/, Si + s/ + - + 1, s/ - Si + - + 1; 2s/ + 1, s/ - Si + - + 1 - n^; 1 , 



(52) 



where 3F2 denotes a generalized hypergeometric function [39]. The generalized hypergeo- 
metric function 3F2 is calculated numerically through a summation of its series represen- 
tation for which a Fortran code is available [ID] . It can be easily checked that if the first 
argument of 3F2 in fl52l) is zero, Uf = 0, one has 3F2 = 1. In this case Eq. fl52|) reduces to 
the following simpler expression 

o(^ r n T^ ^ r(s, + Sf + { + l)r(n, + Si-sj-^) 

bi[ni,li\u,lf) = —rjTT — = ;— . (53) 

^ ^^ c'/2 V^.!r(2s; + l)T{n, + 2s, + l)r(s, - s/ - |) ^ ^ 

This corresponds to a transition from a non-yrast to an yrast state. The integrals for 
the yrast intraband transitions, Eqs. (50) and (51) in [18], are directly obtained from 

19 



Eq. (!53|) when rii = 0. Simple explicit forms of the Si integrals for interband and intraband 
transitions in the particular cases up to n = 2, which are of practical interest, are given 
below 



X 



Siil,h;l,If) 

^ [{2s, + l){2sf + 1) - (s, + Sf--){s, + Sf + - + 1) 
r(s, + s/ + | + l) 



v/r(2si + 2)r(2s/ + 2)' 



(54) 



5/(2,/.;!,// 
V2 
2c 



^^2{s, + l){2s. + l){2sf + l) 



{S^ + Sf + - + 1] 



2{s, + l)(2s, + Asf + 3)-{si + Sf + - + 2){3si + Sf-- + 2) 
T{si + Sf + ^ + l) 



y/r{2s, + 3)T{2sf + 2y 



(55) 



Si{2,If,2,If 
1 



2c '/2 



A{si + l){2s., + l){sf + l){2sf + l] 



- (s. + ^/+2 + l) 



16{si + l){sf + l){si + Sf + I) 



{Si + Sf + - + 2)|2(si + l){2si + 1) + 2{sf + l){2sf + 1) + lQ{si + l){sf + 1) 



- (5. + s/+2+3)(3s. + 3s/-- + 4) 



r(g, + g/ + | + l) 
v/r(2si + 3)r(2s/ + 3)' 



(56) 



Appendix C: Explicit form of the integrals over 



The integrals over the angular variable (p, fj4T]) . with the relevant parities ttj and nj can be 
obtained in the following explicit forms. For A = 2 the integral I2 with ki = k2 = k = odd 
(++) or even ( ) is 



\A:+1 



I^±{k) = -Cat + 



IT 



Ak 



4 ^-^ sin(m7r/2) 

1 + - > — ^ '-^ 

TX ^-^ m 



(57) 
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where Cat = V°^ ^ /n ,\\i ~ 0.915965594177... is the Catalan constant. In the case of 
^1 7^ ^2, both odd or even, the integral is 



(58) 





I^^ih^h) = 


] 

2\k2 


-h\ 


|fc2 — fcl| — 1 . , , , 

4 v-^ sin(m7r/2) 
1 + - > — ^ -^ 

m=l 






+ 


2{k2 + h) 


. k2+ki — l . , ,„N 

4 ■^^ sin(m7r/2) 
TT ^-^ m 

m=l 




A = 


= 3 one has 






l^ {h,h) = 


1 


■(_l)(fe2-fei-l)/2 f_l\ik2+ki~l)/2- 



(59) 



where /ci = 1, 3, 5, . . . , /c2 = 2, 4, 6, . . . For A = 1 the integral is obtained in the form of 
an infinite, but reasonably converging series 



±oo itoo 



n 



O-TT / ^ / ^ / ^ 



sign(— nj 



27r 



m,=±l n=±l u=±l 



\mn\ 



:i-<5, 



sin[(A;2 + z/fci + m + n)|] vr 



k2+uki ,—in—n ) 



{k2 + Z^fci + 772 + n) 



"*" D l^2+'^ki,—m~n 



where ki = 1,3,5, .. . , k2 = 2,4,6, . . 



(60) 
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Table 1: Parameters of the model fits. 



Nucl 


u [MeV/n] 


h [h-^] 


do [h^] 


c 


P 


^Iff N 


i52Sm 


0.295 


2.450 


78.8 


113.2 


0.854 


1.01 


i54Sm 


0.205 


4.625 


108.5 


132.6 


0.808 


1.017 


i54Gd 


0.306 


2.948 


114.7 


113.4 


0.777 


1.048 


i56Gd 


0.439 


1.642 


197.6 


141.5 


0.849 


0.723 


i58Gd 


0.168 


3.626 


42.6 


39.7 


0.864 


0.435 


236U 


0.402 


1.404 


539.3 


343.4 


0.949 


0.134 


looMo 


0.318 


2.674 


1.366 


54.6 


0.715 


0.282 



Table 2: Resulting mass parameters B2 and B^ (in /i^/MeV), Eq. ([3]), and parameters of 
the model potential C2 and C3 (in MeV), (^2 and d^ (in /i^/MeV), Eq. (|1]). The semi-axes 
(sa) of the ellipsoidal potential bottom (3^ and P^, Eq. fH^ . at angular momentum J = 
are given in columns 8 and 9. 



Nucl 


^2 


53 


C2 


c. 


d2 


d^ 


/3f 


PT 


i52Sm 


525 


241 


45.8 


21.0 


429 


197 


0.252 


0.371 


is^Sm 


987 


303 


41.7 


12.8 


427 


131 


0.279 


0.504 


i54Gd 


613 


127 


57.6 


11.9 


416 


86 


0.263 


0.578 


is^Gd 


447 


197 


86.2 


38.0 


545 


240 


0.255 


0.384 


i58Gd 


317 


156 


8.9 


4.4 


175 


86 


0.407 


0.579 


236U 


948 


760 


153 


123 


1351 


1083 


0.226 


0.252 


looMo 


337 


7 


34.0 


0.7 


252 


6 


0.112 


0.759 
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Table 3: Theoretical and experimental values of B(E1), B(E2) and B(E3) transition proba- 
bilities in Weisskopf units (W.u.) for alternating-parity spectra of several even-even nuclei. 
Notations: g (ground-state band), 61 (first /3-band), 62 (second /3-band), nl (first negative- 
parity band), n2 (second negative-parity band), n3 (third negative-parity band). The data 
are taken from [55] except for those for B(i?3; 3~^ — )■ 0^) transitions, which are taken from 
[5B] . The parity signs (-I-) for the even and (— ) for the odd angular momenta, respectively 
are omitted in the labels of the states to avoid overloading of notations. The uncertainties 
(in parentheses) refer to the last significant digits in the experimental data. 



Muh Transition Th [W.u.] Exp [W.u.] Muh Transition Th [W.u.] Exp [W.u. 



152 



Sm 



E2 


2, ^Og 


141 


144 (3) 


E2 


3n2 — ^ ln2 


52 




E2 


4,^2, 


210 


209 (3) 


E2 


5n2 -^ 3„2 


63 




E2 


^^"^9 


248 


245 (5) 


E3 


3nl — ^ 0, 


14 


14(2) 


E2 


Sg^Qg 


284 


285 (14) 


E3 


3„,2 -^ Ofei 


10 




E2 


10, ^ 8, 


322 


320 (3) 


E3 


Inl -^ 4, 


69 




E2 


12, -> 10, 


363 




E3 


ln2 -> 4fei 


70 




El 


Ini -^ 0, 


0.0041 


0.0042 (4) 


E2 


2bl -> 0, 


1.26 


0.92 (8) 


El 


l„i — )■ 2, 


0.0088 


0.0077 (7) 


E2 


4w ^ 2, 


0.2 


0.7 (2) 


El 


3nl — ^ 2, 


0.0056 


0.0081 (16) 


E2 


2fei — )■ 2, 


4.6 


5.5 (5) 


El 


3nl -^ 4:g 


0.0087 


0.0082 (16) 


E2 


4fei -^ 4, 


4.2 


5.4 (13) 


El 


ln2 -^ Ofei 


0.0041 




E2 


2w ^ 4, 


27.4 


19.2 (18) 


El 


1„2 — ^ 2bi 


0.0095 




E2 


461 -^ 6, 


35 


4(2) 


E2 


2w ^ Ofei 


160 


107 (27) 


E2 


Ofei ^ 2, 


30 




E2 


4w — ^ 2fci 


232 


204 (38) 


El 


Inl — ^ 2f,i 


0.00402 


0.00013 (4) 


E2 


3nl — ^ Inl 


47 




El 


Inl — ^ Obi 


0.0023 




E2 


Onl ^ 3nl 


58 


154 


El 
Sm 


ln2 ^ 0, 


0.00006 




E2 


2,^0, 


168 


176 (1) 


E2 


5n2 — ^ 3„2 


82 




E2 


4,^2, 


247 


245 (6) 


E2 


3n3 ~^ ln3 


72 




E2 


6.^4, 


287 


289 (8) 


E3 


3nl -> 0, 


10 


10(2) 


E2 


8,^6, 


322 


319 (17) 


E3 


Inl -> 4, 


50 




E2 


10, ^ 8, 


358 


314 (16) 


E3 


3n2 — ^ Ofci 


77 




E2 


12, ^ 10, 


398 


282 (19) 


E3 


ln2 -^ 4w 


381 




El 


Inl ^ 0, 


0.0051 


0.0058 (4) 


E3 


3n3 — ^ 0b2 


6 




El 


Inl — ^ 2, 


0.0110 


0.0113 (7) 


E3 


InS — ^ 4b2 


62 

continues 


on next page 
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Table 3, continued 



Mult Transition Th [W.u.] Exp [W.u.] Mult Transition Th [W.u.] Exp [W.u. 



El 
El 
El 
El 
El 
El 
E2 
E2 
E2 
E2 
E2 
E2 
E2 



3„i -^ % 

ln2 -^ Obi 
ln2 -^ %l 
ln3 — > 0b2 
ln3 '^ 252 
261 — ^ O51 

4w -> 2. 
2^2 ^ 0, 



■61 



'62 



462 —^262 

<Jnl ^ J-„i 

t'nl ^ <Jnl 

3n2 ~^ ln2 



0.0069 
0.0106 
0.0109 
0.0231 
0.0044 
0.0109 

65 

93 

68 

97 

60 

72 

69 



0.0080 (11) 
0.0092 (13) 



E2 
E2 
E2 
E2 
E2 
E2 
El 
El 
El 
El 
El 
E3 
E3 



061 -> 2, 
261 ^ 0, 
261 -^ 2g 
261 -^ 4^ 

062 -^ 2g 

062 — > 251 
061 -^ Inl 
ln2 ~^ Og 
l?i2 — > 062 
l?i3 — ^ 061 
InS — ^ Og 
3n2 — ^ Og 
3n3 — ^ Og 



1 

0.36 

0.39 

0.27 

5 X 10-^ 

16 

0.0005 

0.0005 

0.0058 

3 X 10-^ 

8 X 10-5 

1.7 

0.4 



12(3) 
<0.58 
<1.3 
<2.4 



154 



Gd 



E2 


2,^0, 


160 


157 


(1) 


E2 


5n2 ~^ 3„2 


64 




E2 


4.^2, 


235 


245 


(9) 


E2 


3?i3 "-^ l?i3 


51 




E2 


6.^4, 


273 


285 


(15) 


E3 


3nl ^ 0, 


21 


21(5) 


E2 


8,^6, 


306 


312 


(17) 


E3 


3n2 — > 061 


32 




E2 


lOg -^ 8g 


340 


360 


(4) 


E3 


3n3 -^ 062 


144 




E2 


12, ^ 10, 


377 






E3 


Inl — ^ 4, 


102 




El 


Inl -> 0, 


0.0102 


0.0436 


E3 


ln2 ^ 4fei 


179 




El 


Inl — ^ 2, 


0.0216 


0.0485 


E3 


ln3 -^ 4f,2 


708 




El 


3„i — )■ 2, 


0.0137 






E2 


061 ^ 2, 


25 


52 (8) 


El 


3nl — > 4, 


0.0207 






E2 


261 ^ 0, 


1.23 


0.86 (7) 


El 


ln2 -> 061 


0.0152 






E2 


261 ^ 4, 


22.6 


19.6 (16) 


El 


l?i2 — ^ 261 


0.0333 






E2 


062 ^ 2, 


0.0553 




El 


ln3 ^ 062 


0.0333 






E2 


062 — ^ 2^,1 


14 




El 


ln3 — > 262 


0.0706 






El 


Inl ^ 061 


0.0054 


0.0057 


E2 


261 ^ Ofei 


177 


97 1 


(10) 


El 


Inl — > 261 


0.0099 


0.0064 


E2 


461 -^ 261 


256 






El 


ln2 -^ O3 


2x10-5 




E2 


262 -> 062 


85 






El 


ln2 ^ 062 


0.0094 




E2 


462 -> 2b2 


122 






El 


ln3 -^ 061 


0.00023 




E2 


3nl — )■ Inl 


55 






El 


ln3 ^ O3 


2x10-6 




E2 


5nl — )■ 3„i 


67 






E3 


3n2 — ^ O3 


1.8 




E2 


3n2 — ^ ln2 


54 






E3 


3n3 — ^ Og 


0.05 
continues 


on next page 
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Table 3, continued 



Mult Transition Th [W.u.] Exp [W.u.] Mult Transition Th [W.u.] Exp [W.u. 



156 



Gd 



E2 


2g^0g 


150 


187 (5) 


E2 


E2 


4,^2, 


219 


263 (5) 


E2 


E2 


69^4, 


249 


295 (8) 


E3 


E2 


8g^6g 


273 


320 (14) 


E3 


E2 


10, ^ 8, 


296 


314 (14) 


E3 


E2 


12, -> 10, 


321 


300 (3) 


E3 


El 


Inl -^ 0, 


0.0006 


0.0019 (14) 


E2 


El 


Inl — ^ 2, 


0.0013 


0.0025 (18) 


E2 


El 


3nl — )■ 2, 


0.00083 


0.00098 (21) 


E2 


El 


3nl — ^ 4, 


0.0012 


0.00077 (16) 


E2 


El 


ln2 — ^ Obi 


0.0013 




E2 


El 


1„2 — ^ 2bi 


0.0026 


0.0005 (3) 


El 


El 


3„2 — > 2bi 


0.0016 




El 


E2 


2w ^ Ow 


74 


52 (23) 


El 


E2 


4bi — )■ 2fei 


107 


280 (15) 


El 


E2 


651 ^ 461 


120 




El 


E2 


3nl — ^ Inl 


46 




E3 


E2 


5nl — > 3„i 


56 







3n2 -^ 1 



-n2 



5n2 —^ 3„2 
3ral — ;■ 0, 

3n2 — ^ Obi 
Inl -^ 4, 

ln2 ^ 4, 

Obi -^ 2, 
261 -> 0, 



^61 



'61 



461 ^ 2, 
461 -> 6, 
261 ^ 4, 

Inl -^ Of 
ln2 ^ 0, 
ln2 — ^ 2, 
3n2 — ^ 2, 
3n2 ^ 4, 
3n2 ~^ 0, 



44 
53 
16.9 
64 
73 
282 

5 

0.32 

0.1 

5.6 

4.3 

0.0002 

6 X 10*6 

2 X 10-5 

5 X 10-6 

2 X 10-5 

0.21 



16.9 (7) 



8(4) 

0.63 (6) 

1.3 (7) 

2.1 (11) 

4.1 (4) 

0.0004 (3) 

0.0019 (7) 

0.0043 (15) 

0.0019 (14) 

0.0031 (4) 



158 



Gd 



E2 


2,^0, 


181 


198 (6) 


E2 


E2 


4,^2, 


274 


289 (5) 


E2 


E2 


63^4, 


332 




E2 


E2 


8,^6, 


393 


330 (3) 


E2 


E2 


10, ^ 8, 


460 


340 (3) 


E2 


E2 


12, ^ 10, 


532 


310 (3) 


E2 


El 


Inl -> 0, 


0.0001 


9.8443x10-5(4) 


E2 


El 


Inl — ^ 2, 


2.5x10-4 


9.6515x10-5(6) 


E2 


El 


3nl -^ 2, 


0.00015 


0.00033 (10) 


E2 


El 


3nl — ^ 4, 


0.00028 


0.00029 (8) 


E2 


El 


5nl — ^ 4, 


2.02x10-4 


7.4324x10-4(13) 


E2 


El 


5nl -^ 6, 


3.62x10-4 


5.8691x10-4(8) 


El 


El 


3„2 — > 261 


0.00011 


> 0.00035 


El 


El 


ln2 — > 061 


8.02x10-5 




El 


El 


ln2 — ;■ 2bi 


0.0002 




El 



061 -> 2, 
261 ^ 0, 
261 -^ 2, 
461 ^ 4, 

261 -^ 4, 

062 ^ 2, 

262 -> 0, 
262 -^ 4, 
461 -> 2, 
461 -> 6, 

062 — > 2bi 
Obi — ^ Inl 
2bl — ^ Inl 
2bl — ^ 3„i 
ln2 -^ 2, 



1.1652 

0.31 (4) 

0.079 (14) 

0.37 
1.39 (15) 

2.09 
0.37 (4) 
0.38 (6) 

1.32 

3.16 

3.314x10-6 

6.4x10-5(8) 

1.89x10-4(24) 

0.0064 

continues on next page 



8.7619 

2.36 

2.913 

2.40 

2.96 

1.86 

0.68 

0.43 

3.75 

1.30 

57 

2.7x10 



-5 



8.3x10-6 
2x10-5 
4x10-5 
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Table 3, continued 



Mult Transition Th [W.u.] Exp [W.u.] Mult Transition Th [W.u.] Exp [W.u.] 



El 


ln3 -^ 


062 


0.0004 




El 


El 


ln3 -^ 


2^2 


0.0009 




El 


El 


3„3-> 


2^2 


0.0005 




El 


E2 


2w^ 


061 


200 




El 


E2 


4w -^ 


261 


288 


455 


El 


E2 


2b2^ 


062 


217 




El 


E2 


4^2 -^ 


262 


308 




El 


E2 


3„i^ 


1„1 


185 




El 


E2 


5„i -^ 


3nl 


227 


369 (6) 


El 


E2 


3„2^ 


ln2 


200 


> 1600 


El 


E2 


5„2 -^ 


3n2 


240 




El 


E2 


3„3^ 


ln3 


241 




El 


E3 


3nl — >■ 


o« 


11.9 


11.9 (7) 


El 


E3 


Inl ^ 


4« 


81 




El 


E3 


3„2^ 


061 


519 




E3 


E3 


3„3 ^ 


06. 


102 




E3 



ln2 -> Og 
3n2 ~-^ 2g 
3n2 -^ 4g 
062 ->1 



nl 



262 -^ Inl 

262 — ^ 3ni 

062 ^ 1 

262 -> 1 

262 — > 3, 



n2 
n2 



'n2 



461 -> 3„i 

461 -^ 5ni 
ln2 — ^ 062 
ln3 ""^ 061 
ln3 — ^ Og 
3n2 — ^ O3 
3n3 — ^ Og 



2x10-5 
3x10-5 
3x10-5 

2x10-^ 

2x10-^ 

2x10"^ 

6x10-5 

1.8x10-5 

4.2x10-5 

7.7x10-6 

2.1x10-5 

2x10-5 

3x10-6 

0.00001 

5 

2 



0.0035 (12) 

>0.0011 

>0.0015 

5.7831x10-5 

2.7x10-6(19) 

3.7x10-5(5) 

6.02x10"^ 

1.50x10-^(21) 

2.40x10-^(5) 

4.63x10-^ 

6.12x10-^ 



236 



U 



E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
E2 
El 
El 
El 
El 
El 



2,^0, 
4.^2, 
6.^4, 

8g ^ Qg 
lOg -^ 8g 

12g -^ lOg 

Ug -^ 12g 

16g -^ Ug 

18g -^ 16g 

20g -^ 18g 

22g -^ 20g 

2Ag ^ 22^ 
26^ ^ 24^ 

28g -^ 26g 

Inl -^ Og 

J-nl ' ^g 
<Jnl ^ '^g 
3nl — ^ 4:g 
l?i2 — ^ 061 



237 

342 

382 

408 

429 

450 

471 

493 

516 

539 

564 

590 

617 

645 
2.7x10-^ 
5.5x10-^ 
3.5x10-^ 
4.8x10-8 
2.0x10-8 



250 (10) 

357 (23) 

385 (22) 

390 (4) 

360 (4) 

410 (7) 

450 (5) 

380 (4) 

490 (5) 

510 (8) 

520 (12) 

670 (13) 

670 (19) 

1100 (5) 

2.7x10-8(4) 



E2 
E2 
E2 
E2 
E2 
E2 
E2 
E3 
E3 
E3 
E3 
E2 
E2 
E2 
E2 
El 
El 
El 
E3 



261 — ^ 061 

461 ^ 2, 



'61 



3nl y Inl 

5ni y 3„i 

7nl -^ 5 
3n2 -^ 1 
5n2 ~^ 3- 



'nl 
n2 



'n2 



Inl -^ 4:g 
3nl — ^ Og 
ln2 -^ 4fei 
3n2 — > 061 

061 ^ 2g 

261 ^ 0^ 
461 -^ 2^ 
261 -^ 4^ 

061 — > Inl 
261 — > Inl 
ln2 — ^ Og 
362 ^ 0, 



112 

160 

68 

80 

87 

54 

64 

62 62 (9) 

15 23 (3) 

695 

172 

6 

0.66 

0.59 

4 

1.2x10-8 

4.6x10-9 

1.6x10-9 

0.14 

continues on next page 
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Table 3, continued 



Mult Transition Th [W.u.] Exp [W.u.] Mult Transition Th [W.u.] Exp [W.u. 



El 1„9^2m 4.0x10-^ 



100 



Mo 

E2 2g -^ Og 22.7 37.0 (7) E2 3„i -^ l„i 16 

E2 Ag -^ 2g 50 69 (4) E2 5„i -^ 3„i 21 

E2 6g -^ Ag 84 94 (14) E2 7„i -^ 5„i 26 

E2 8g -^ 6g 120 123 (18) E2 3„2 -> 1„2 18 

E2 IO3 -^ 8g 156 E2 5„2 -^ 3„2 22 

El Ini ^ O3 2x10-6 E3 3„i ^ O3 34 34(3) 

El lni^2g 1x10-5 E3 3„2 ^ Obi 5 

El 3ni^2g 7x10-6 2.7x10-6(9) E3 l„i ^ 4^ 899 

M 3„i^43 2x10-5 E2 0bi-^2g 72 92(4) 

El 1„2 ^ Ow 2x10-^ E2 2bi^0g 0.5 0.62(5) 

El 1„2 ^ 2w 1x10-5 E2 4:bi^2g 3 

El 3„2 ^ 2w 3x10-6 El l„i ^ Obi 8x10-6 

El 3„2 ^ 4w 3x10-5 El l„i ^ 2bi 2x10-5 

E2 2bi ^ Ofei 25.4 5.5(8) El 3„i ^ 2^1 1.4x10-5 2.5xl0-5(^ 

E2 4bi ^ 2f,i 45 El 1„2 ^ 0^ 5x10-^ 

E2 6bi -^ 4fei 75 E3 3„2 ^ Og 22 
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Figure 1: (Color online) Theoretical and experimental alternating-parity bands in ^^^Sni. 
Data from |3l]. The oscillation quantum numbers n, kn and kn are given above the 
theoretical bands. 
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Figure 2: (Color online) The same as in Fig. 1, but for Sm 
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Figure 3: (Color online) The same as in Fig. 1, but for Gd. 
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Figure 4: (Color online) The same as in Fig. 1, but for Gd. 



33 



2000 



1500 



> lOOOi 
a> 

Oi 

c 

LU 



n=2 



158 



Gd 



n=0 



-1732.53 n=1 



^.1^5161 5-1^63^0^^1540.19 
-1457.07 Q+lJ26s2p 



J- 1995.91 
,-1905.89 



+ 168332 



10 



4278.50 



^+ 1347.54 - 



- 1393.69 
-1294.16 



g- 1228.12 ^ 1204.4J 







+ 1140^2 



8 



+ 912.76 



-1057.91 



958.38 



-1978.04 



-1856.32 



1684. 



,+ 1636. - 1639.34 4^^-^^^^ 



10 



+ 1350. 



_- 1391. 4+ 1406.7 J- 1402.94 ' 

^+ 1259.87 1- 1263.51 
,-1176.48 n+U26:l? 



.+ 1517.48 
.+ 1452.35 



- 1041.64 
-977.145 



8 



+ 904.12 



500 



.+ 580.30 



^+ 539.021 



^+ 299.07 



.+261.457 



.+ 95.41 



,+79.5128 



0+' 



theory 



0+' 



experiment 



Figure 5: (Color online) The same as in Fig. 1, but for Gd. 
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Figure 6: (Color online) The same as in Fig. 1, but for U. 
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Figure 7: (Color online) The same as in Fig. 1, but for Mo 
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Figure 8: (Color online) Density distribution Pnki{.P2, Pz) = l'^nfc/(/^2? /^s)!^ for: (a) A; = 1, 
1 = 2 and (h) k = 2, I = 1 at n = with schematic parameters (see the text). The model 
space corresponds to the (32 > half-plane. 



37 



«a 



-0.2 



qs; 0.0 




<=5. 



qS 0.0 - 



-0.4 -0.2 0.0 0.2 0.4 




Figure 9: (Color online) Contour plots of the density distribution Pnki{.P2i /^s) for: (a) A; = 1, 
/ = 2, (b) A; = 2, / = 1, (c) fc = 3, / = 2 and (d) fc = 4, J = 1 at n = with the schematic 
parameters (see the text). The ellipsoidal curves outline the potential bottom. The model 
space corresponds to the /32 > half-plane. 
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